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ABSTRACT 

The WMAP satellite has provided high resolution, high signal to noise ratio maps of 
the sky in five main frequency bands ranging from 23 to 94 GHz. These maps consist 
in noisy observations a mixture of Cosmic Microwave Background (CMB) anisotropics 
and of other astrophysical foreground emissions. We present a new foreground-cleaned 
CMB map, as well as a new estimation of the angular power spectrum of CMB tem- 
perature anisotropies, based on 7 years of observations of the sky by WMAP. The 
method used to extract the CMB signal is based on an implementation of minimum 
variance linear combination of WMAP channels and of external full-sky foreground 
maps, on a frame of spherical wavelets called needlets. The use of spherical needlets 
makes possible localised filtering both in pixel space and harmonic space, so that the 
ILC weights are adjusted as a function of location on the sky and of angular scale. 
Our CMB power spectrum estimate is computed using cross power spectra between 
CMB maps obtained from different individual years of observation. The CMB power 
spectrum is corrected for low-level biases originating from the ILC method and from 
foreground residual emissions, by making use of realistic simulations of the whole 
analysis pipeline. Our error bars, compatible with those obtained by the WMAP col- 
laboration, are obtained from the combination of two terms: the internal scatter of 
individual Ci in each £ bin, and a term originating from uncertainties in our correction 
for biases due to empirical correlations between CMB and foregrounds, as well as to 
residual foregrounds in the CMB maps. Our power spectrum is essentially compatible, 
within error bars, with the result obtained by the WMAP collaboration, although it 
is systematically lower at the lowest multipoles, more than expected considering that 
the two estimates are based on the same original data. Exhaustive investigations of 
the presence of a possible bias in our estimate fail to explain the difference. Compari- 
son with several other analyses confirm the existence of differences in the large scale 
CMB power, which are significant enough that until the origin of this discrepancy is 
understood, some caution is recommended in scientific work relying much on the exact 
value of the CMB power spectrum in the Sachs- Wolfe plateau. 
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1 INTRODUCTION 

Cosmic Microwave Background (CMB) anisotropies provide 
a snapshot of what the universe looked like at the moment 
of recombination. The statistical properties of CMB fluctu- 
ations depend on the original primordial perturbations from 
which they arose, as well as on the subsequent evolution of 
the universe as a whole. This makes the precise measure- 
ment of CMB power spectra, both in temperature and po- 
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larisation, a gold mine for understanding and describing the 
universe from early times to now. 

For cosmological models in which initial perturbations 
are Gaussian, the information carried by CMB anisotropies 
is completely characterised by the multivariate angular 
power spectrum of the observations (i.e. the auto and cross 
power spectra of temperature and polarisation maps). The 
angular power spectrum of the CMB is sensitive to the val- 
ues of the different cosmological parameters which primar- 
ily describe the fundamental properties of the universe, such 
as its matter content, age, expansion history, global geome- 
try and the properties of the initial fluctuations that seeded 
the large scale structure observable today. The possibility to 
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use this for constraining the cosmological parameters hinges 
upon the existence of acoustic oscillations in the primor- 
dial plasma before last scattering of CMB photons, and on 
the ability of CMB experiments to disentangle such primor- 
dial fluctuations from those induced by cosmic structures at 
lower redshift. This is why a key goal of the CMB commu- 
nity is to measure the true CMB temperature anisotropies 

and polarisation in the sky. 

Since the discovery of the CMB by iPenzias fc WilsonI 

tremendous efforts on the instrumental side have 
been made to improve sensitivity and angular resolution 
of CMB experiments. Recently, multi-frequency data for 
the complete sky have become available fro m 7 years of 
obser vation with the WMAP space mission jjarosik et al.l 
|2010| ). WMAP confirmed the leading theory in cosmology, 
the ACDM (cold dark matter) model of a universe governed 
by Einstein's theory of general relativity, and the evolution 
of which is dominated by the impact of a cosmological con- 
stant A, equivalent to a contribution, in the total energy 
density, of vacuum energy, or more gene rally 'dark energy' 
IILarson et al ] |201ll : lKomatsu et al.|[2O10l ). 

The WMAP data has already allowed a determination 
of several central cosmological parameters with good accu- 
racy. These determinations are based on a number of sim- 
plicity assumptions, which need to be tested by more ac- 
curate and extensive measurements. The Planck satellite of 
the European Space Agency, next generation CMB space 
mission after WMAP, surpasses its pr edecessor in resolu- 
tion, sensitivity and frequency cov erage (jTauber et al.ll2010l : 
iPlanck Collaboration et al.ll20lil ). and will provide an up- 
date of the present picture and of cosmological parameter 
estimates in early 2013. 

Determining the cosmological parameters with CMB 
experiments, however, requires careful cleaning of the CMB 
maps from contamination by galactic and extra-galactic 
foregrounds. In the low frequency microwave regime (below 
about 100 GHz) the strongest contamination comes from 
galactic synchrotron and free-free emission. At higher fre- 
quencies, where synchrotron and free-free emissions are low, 
dust emission dominates. In addition, extragalactic point- 
like sources are a significant contaminant at high galactic 
latitude and on small scale. Even when the brightest sources 
are identified and subtracted-out or masked, their residual 
contribution may dominate over the cosmic variance uncer- 
tainty on small angular scales, an d bias the measurement 
of the CMB power spectrum (see iDelabrouille fc Cardoso! 
l|2009h for a review on component separation in CMB obser- 
vations) . 

Most of the foreground emissions are strongly corre- 
lated between frequency channels. Except for the very faint 
kinetic SZ signal, the underlying frequency dependences of 
the emissions differ from that of the CMB anisotropies. 
Exploiting this facts, a model-independent method to re- 
move foregrounds from the multi-frequency observations of 
CMB, the so-called 'Internal Linear Combination' (ILC), 
has been proposed to extract the CMB signal from the multi- 
frequency data such as that o f WMAP or Planck (see, e.g., 
iTegmark fc Efstathioul l| 19961 )). The idea behind the ILC 
method is to find the linear combination of the available 
maps which has minimal variance while retaining unit re- 
sponse to the CMB. The main advantage of this foreground 
cleaning method is that it does not require any assump- 



tion about the foregrounds. Another advantage is that it 
is easy to implement, and computationally fast. Finally, it 
can be extended to impose the rejection of a particular fore- 
ground if needed (jRemazeilles et al. 2011a,), or to extract the 
total combined emission of several correlated foregrounds 
(|Remazeilles et al.|[2011bl ). 

There are, however, also drawbacks t o the method. 
As d is cussed by a num b er of authors (Hinshaw et al.l 
20071: ISaha et al. I l2008l IDelabrouille fc Cardosd l2009l : 



Delabrouille et all l2009h . the component of interest (the 



CMB in our case) and foreground signals must be uncor- 
related for proper ILC performance. This can be, on finite 
data sets, only approximately true, and empirical correla- 
tions between the CMB and foregrounds genera te a bias in 
the re constructed CMB. In addition, as shown bv lDick et al.l 
l|2010t ). the ILC method tends to amplify calibration errors, 
and reconstruct a CMB map which can be largely under- 
calibrated, in particular for high signal-to-noise ratio obser- 
vations. These sources of error must be monitored carefully 
in the analysis of CMB maps obtained in this way. 

In this paper, we address the problem of measuring 
as precisely as possible the CMB temperature power spec- 
trum from WMAP 7-year observations, with close to min- 
imal noise and contamination by foreground emission. The 
paper is organised as follows: In section [2] we describe the 
methodology to estimate the CMB and its power spectrum 
using an ILC on wavelet decompositions of sky maps ob- 
served at different frequencies. The implementation of this 
on WMAP 7-year data is described in section [3] and section 
|4] The results are discussed in section [5] We conclude in 
section [S] 



2 THE NEEDLET ILC 

Component separation with an ILC method can straight- 
forwardly be performed in real space or in harmonic space. 
By this we mean that different ILC weights can be com- 
puted in different regions of the sky, or in different regions 
of harmonic space. This allows for variations of the data 
covariance matrix in either space. 

The ILC in harmonic space however does not take 
into account the fact that noise is the dominant source of 
CMB measurement error at high galactic latitude while fore- 
ground signals dominate at low galactic latitude. On the 
other hand, the ILC in pixel space does not take into ac- 
count the fact that noise dominates at high angular fre- 
quency (small scales) while foreground emission dominates 
on large scales. In order to overcome this problem, we follow 
the approach of lDelabrouille et al.l l|2009l l and implement the 
ILC on a frame of spherical wavelets, called needlets. This 
special type of spherical wavelets allows localised filtering 
in both pixel space and harmonic space because they have 
compact support in the harmonic d omain, while still being 
very well localised in pixel domain ([Narcowich et al.l [20061 : 
iMarinucci et al.l l2008l : ICuilloux et all 120091 '). Needlets have 
already been used in various analyses of WMAP data be- 
sides compone nt separation and power spectrum estimation, 
for instance bv lPietrobon et al] (120081 to detect features in 
the CMB, and bv lRudiord et al.l ([2009i' ) to put limits on the 
non-gaussianity parameter /nl- 

Various versions of the ILC on WMAP data have been 
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implemented by a number of authors so far llTegmark et al 



20031 : [Eriksen et al] 120041: iHinshaw et ai] |2007l: iPark et al 
200?; ^Saha ct al."2008': Jpelabrouille et al.l l2009l : lKim et al 

2009; Souradeea .201% 
tion on WMAP data 
endent component a nalysis 



The filters hj is chosen in such a way that 



Alternate component separa- 
has been performed using inde- 
methods such as SMICA 



= 1, 



(7) 



which permits us to directly reconstruct the original maps 
from their filtered maps a"^^ using the same set of filters. 



Patanchon et al.| |2005|) C^A jBonaldi et al |2007|). and ^erms of hj , the spherical needlets are defined as. 



FASTICA (Main o et al] 120071 : iBottino et £01 12008'. 20fO^ 



Specific analyses to extract foreground component emissions 
have permitted to isola te the total foreground contamina- 
tion (Chosh et al.ll201ll) or detect the thermal SZ emission 
of ga laxy clusters ( Atrio-Barandela et al.|[200^ : lMelin et al.l 



2.1 The data model 

A CMB space mission such as WMAP provides full- 
coverage, multi-frequency anisotropy maps of the sky 
jtObs,c^^-j different frequency bands (channels). The 

observed signal T'~'^^''^ (n) in channel c can be modelled as, 



nOBS,( 



'{h) = / dnf,' B"(n,n') T^'°'"(n')-f rN'"(n), (1) 

where T^^'^''^{h) is the signal (sky) component, itself decom- 
posed in the sum of CMB and foreground components, 

T ' (n) = a T (n)+T (n), (2) 

a'^ being the CMB calibration coefficient for the channel c. 
Up to calibration uncertainties, a° = 1 for all WMAP chan- 
nels. If, in addition to WMAP data, we use extra ancillary 
data which serve as foreground templates to help foreground 
subtraction, as done in the present work, the coefficients a" 
vanish for such data sets. 

The beam function B'^{n,h') represents the smoothing 
of the signal due to the finite resolution of the observations. 
The beam is assumed here to be circularly symmetric, i.e., 
B'^(n, n') depends only on the angle 9 = cos~^(n.n') between 
the directions h and fi . We may then expand this function 
in terms of Legendre polynomials. 



B'^{h, n) 



^ Bi Pi{n.n 



47r 



(3) 



The term T^''^(n) in equation[T]represents the detector noise 
in channel c. Unlike the CMB and foreground components, 
instrumental noise is not affected by the beam function. 

Equation [1] can be recast, in the spherical harmonic rep- 
resentation, as: 



OBS.c c DC CMB , DC FG,c , N 



(4) 



(8) 



where {£,jk} denote a set of cubature points on the sphere, 
corresponding to a given scale j. In practice, we identify 
these points wit h the pixel centres in the HEALPix pix- 
elisation scheme (|G6rski et al.|[2005l ). The cubature weights 
\jk are inversely proportional to the number Nj of pixels 
used for the needlet decomposition at scale j, i.e. Ajfe = 4r"- 
The needlet coefficients for CMB temperature anisotropies 
T{n) = Y!m=^i o-imYimin) are denoted as, 



ft-fc = / T(n)*jfc(n)df7ft 

!max I 

= hj Biai^ Yimi^jk)- 

1 = m=-l 



(9) 



The linearity of the needlet decomposition implies that the 
needlet coefficients Pjf. corresponding to the filtered map ob- 
tained from the harmonic coefficients a"^^ are a linear com- 
bination of the needlet coefficients of individual components 
and noise at HEALPix grid points ^jk'- 



pC C qCMB . qFG,C . qN,C 



(10) 



where, 



oCMB 

Pjk 



h\Biai^ Yimi^jk) 



1 = m = -l 
'max I 



^.^E E hlB.a^.^'^YUi^k) 



1 = m = ~l 
'max I 
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1 = m = -l 



(11) 



2.3 Implementation of the needlet ILC 

The ILC estimate of needlet coefficients of the cleaned map 
is obtained as a linearly weighted sum of the needlet coeffi- 
cients 



2.2 Implementation of the needlet transform 

Considering that each channel observes the sky at a differ- 
ent resolution, the maps are first convolved /deconvolved, in 
harmonic space, to the same resolution: 

c Bl OBS.C ,r\ 

aim = «im ■ (5) 

Each of these maps af^ is then decomposed into a set of 
filtered maps aj^^ represented by the spherical harmonic co- 
efficients. 



(6) 



oNILC _ c oc 

Pjk - ^jk Pjk 



(12) 



where lo'^^. is the needlet weight for the scale j and the fre- 
quency channel c at the pixel k. Under the assumption of 
de-correlation between CMB and foregrounds, and between 
CMB and noise, the empirical variance of the error is min- 
imum when the empirical variance of the ILC map itself 
is minimum. The condition for preserving the CMB signal 
during the cleaning is encoded as the constraint: 



E' 



^jk 



1. 



(13) 
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The resulting needlet ILC weights Qj^. that minimise the 
variance of the reconstructed CMB, subject to the constraint 
that the CMB is preserved, are expressed as: 



R 



the clei 



(14) 



The NILC estimate of the cleaned CMB needlet coefficients 
is: 



oNILC _ aCMB . 
Pjk — Pjk + 



where _RIJ 



ance matrix iij^ = 



,jf. are empirical estimates of the elements of covari- 
PjkPjk^ for scale j at pixel k. Those 
estimates are obtained each as an average of the product of 
the relevant computed needlet coefficients over some space 
domain Dk centred at k. In practice, they are computed as 



R-Tk = — 5Z ™^ '^')l^lkl3]k , 



(15) 



where the weights Wj{k,k') define the domain "Dk. A sensi- 
ble choice is for instance Wj{k, k') = 1 for k' closer to k than 
some limit angle, and Wj{k,k') = elsewhere, or alterna- 
tively, Wj(k,k') shaped as a Gaussian beam of some given 
size that depends on the scale j (which is what we do here). 

Finally, the NILC estimate of the cleaned CMB map 
can be reconstructed from cleaned CMB needlet coefficients 
using the same set of filters that was used to decompose the 
original maps into their needlet coefficients. The NILC CMB 
temperature map is then 



ENILC -.^ I 
aim Yim(n) 



E( n C 



RFG 
aim 



m 

' Yim{n), (16) 



where the harmonic coefficients residual foreground (af^^) 



and residual noise(a^) 



are given by: 



and 



Oln 



Jk) 



(17) 



(18) 



Equation [16] implies that the NILC estimate of CMB con- 
tains some residual foreground and noise contamination. 



3 WMAP 7 YEAR NEEDLET ILC MAP 

The WMAP satellite has observed the sky in five frequency 
bands denoted K, Ka, Q, V and W, centred on the frequen- 
cies of 23, 33, 41, 61 and 94 GHz respectively. After 7 years of 
observation, the released data includes maps obtained with 
ten difference assemblies, for 7 individual years. One map is 
available, per year, for each of the K and Ka bands, two for 
the Q band, two for the V band and four for the W band. 
These sky maps are sampled using the HEALPix pixelisa- 
tion scheme at a resolution level (nside 1024), corresponding 
to approximately 12 million sky pixels. 



Table 1. List of needlet bands used in the present analysis. 



Band index 




^peak 


^max 


nside 


1 








50 


32 


2 





50 


100 


64 


3 


50 


100 


150 


128 


4 


100 


150 


250 


128 


5 


150 


250 


350 


256 


6 


250 


350 


550 


512 


7 


350 


550 


650 


512 


8 


550 


650 


800 


512 


9 


650 


800 


1100 


1024 


10 


800 


1100 


1500 


1024 



In a first step, for each frequency band, we average 
all the difference-assembly maps obtained at the same fre- 
quency. This yields 5 band-averaged maps. However, con- 
trarily to the 7-year average maps at nside = 512 also pro- 
vided by WMAP, these maps at nside = 1024 are not offset 
corrected. In order to determine the values of the offset for a 
particular frequency band, we have used standard resolution 
7-year band-average maps of WMAP as references. Offset 
values, for all frequency bands, are obtained by determin- 
ing the mean of the difference between the band-averaged 
map being considered (degraded to nside = 512), and the 
released 'reference' map. We implement the needlet ILC on 
these offset-corrected maps. 

In addition to these maps, we use three foreground tem- 
pla tes in our ana lysis: dust at 100 microns, as obtained 
b y ISchlegel et al.l ll QQS). the 408 MHz synchrotron map 
of iHaslam et al.l ~('T981). and the composite all-sky H-alpha 
map of iFinkbe iner (200^ . The offset-corrected WMAP 
maps and fore-mentioned foreground templates are con- 
volved/deconvolved, in harmonic space, to a common beam 
resolution (corresponding to that of the W frequency chan- 
nel of WMAP 7-year release) before the implementation of 
the needlet ILC on the map set. 

Each of these maps is decomposed into a set of needlet 
coefficients. For each scale j, needlet coefficients of a given 
map are stored in the format of a single HEALPix map at 



degraded resolution. The filters ft; used to compute filtered 
maps are shaped as follows: 



for li 



^ peak' 



for I — Ipeak ^ 



for <r I < 

'^"^ ''peak ^ ' ^ f-max 



For each scale j, the filter has compact support between the 
multipoles Z^-^j and ll^ax with a peak at V^^^^f, (see figure 
[1] and table [l}. The needlet coefficients /3jfc are computed 
from these filtered maps on HEALPix grid points with 
resolution parameter nside equal to the smallest power of 2 
larger than Z^^^/S. 

The estimates of needlet coefficients covariance matri- 
ces, for each scale j, are computed by smoothing all possi- 
ble products of needlet coefficient maps PjkPjk with Gaus- 



A needlet ILC analysis of WMAP 7-year temperature data 5 




Multipole{l) 



Figure 1. Needlet bands used in the present analysis. The black 
line shows the normalisation of the needlet bands, i.e. the total 
filter applied to the original map after needlet decomposition and 
synthesis of the output map from needlet coefficients. 




-0.50 ^^^^m 5U mK_CMQ 



4 CMB ANGULAR POWER SPECTRUM 

Component separation with a needlet ILC is the first step 
in our analysis. The cleaned CMB map, however, still con- 
tains residuals of noise and foreground emission, and is im- 
pacted by the ILC bias due to empirical correlations of CMB 
and foregrounds as well as by the effect of calibration un- 
certainties and beam uncertainties. Residual contamination 
is clearly visible, for instance, along a narrow strip on the 
Galactic plane of the map displayed in figure [S] 

In the following, we discuss these various types of errors, 
and our strategy to minimise or evaluate their impact on the 
estimation of the CMB power spectrum. 



4.1 Impact of instrumental noise 

It would be possible to use the estimated noise level in the 
WMAP maps to de-bias the power spectrum of the map 
obtained in section O from instrumental noise contribution. 
This, however, requires, in particular at high multipoles, a 
very accurate estimate of the noise level. 

Thus, instead of trying to de-bias from instrumental 
noise a single power spectrum (computed on a single NILC 
map averaging all seven years of data), we produce inde- 
pendent CMB maps for each of the individual 7 years of 
observations, and compute the CMB power spectrum as an 
average of the 21 cross power spectra obtained from the 7 
clean CMB maps. In practice, all 7 maps are obtained us- 
ing the same set of needlet weights, determined using the 
co-added 7 year observations. Instrumental noise, being un- 
correlated from year to year, does not bias the average power 
spectrum, and can be ignored in the estimate. 



Figure 2. The NILC estimate of CMB temperature anisotropics 
obtained by implementing the NILC on the 7-year band average 
maps (at nside = 1024) provided by the WMAP collaboration. 



sian beams. In this way, an estimate of needlet covariances 
at each point k is obtained as a local, weighted average of 
needlet coefficient products. The process is summarised by 
the diagram: 

{l3jkPjk}k^l...,Npi, - > {"i^m }i = 0,4„^,m=-;,..,0,..,I 
{-Rjfc }fc=l...,JVp>, < i^i°''^}l^O,li,^^.m^-l...,0,..,l 

The full width at half maximum (FWHM) of each of the 
Gaussian windows used for this purpose is chosen to ensure 
the computation of the statistics by averaging about 1200 
samples or more. Choosing a smaller FWHM results in ex- 
cessive error in the covariance estimates, and hence excessive 
bias. Choosing a larger FWHM results in less localisation, 
and hence some loss of effectiveness of the needlet approach. 

Using these covariance matrices, the ILC solution for 
each scale is implemented locally to get the needlet weights 
and hence the estimated CMB needlet coefficients. Finally, 
a full sky CMB map, displayed in figure [2l is synthesised 
from these estimated needlet coefficients. 



4.2 Impact of residual foregrounds 

As the ILC weights used for all years are the same, residual 
foreground emission will be the same in all maps. These 
residuals, although small (as compared to original fore- 
ground contamination in the maps of the individual differ- 
encing assemblies), are thus 100% correlated between the 
different single year maps, and bias the CMB power spec- 
trum. They are dealt with using a combination of masking 
and de-basing, as follows. 

First of all, the impact of most of the contamination 
from residuals of bright point sources is removed from our 
estimate of power spectra by applying to our output CMB 
maps, for each individual year, the point source mask pro- 
vided by the WMAP collaboration, and fiUing-in the holes 
by an interpolation procedure using the values of CMB 
anisotropics in the neighbouring unmasked pixels. We start 
with the masked border pixels, assign to each of them the 
average of all observed pixels which are within a distance of 
two pixel sizes. Then we iterate, increasing the distance for 
averaged pixels by two pixel sizes at each iteration. While 
not theoretically optimal, this procedure works very well in 
practice. 

The impact of galactic contamination on our power 
spectrum estimate is limited by the use of a conservative 
galactic mask, Mg^ (applied after the NILC, so that we 
still get full sky maps). Our mask simply excludes sky re- 
gions in which galactic residuals after component separation 
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6 (in degrees) 

Figure 3. Apodisation kernel for Galactic mask used in our anal- 
ysis 



can be too strong for precise CMB power spectrum estima- 
tion. The apodised mask used in practice is: 

for 0° < < 9r, 
{^)^] for s; |e| < + 

1 for \e\ > Or + eu, 



Ms B 



where 6 = (f ~ ^) is the latitude. The apodisation of 
the mask borders hmits the ahasing of large scales into small 
scales. 



4.3 Validation on simulations and de-biasing 

The method is tested, tuned, and validated using numer- 
ical simulations based on the Planck Sky Model (PSM). 
A general description of the PSM can be found on the 
P SM web pagen as well as in iLeach et al.l (|2008l ) and/or 



in iBetoule et al.l (|2009l ') . The simulated sky emission com- 
prises: 

• a Gaussian CMB temperature map, with a power spec- 
trum matching the WMAP best fit; 

• galactic emission, consisting in the superposition of syn- 
chrotron, free-free, thermal dust, and spinning dust; 

• SZ effects, both thermal and kinetic; 

• a population of point sources. 

Noise maps for the individual years are randomly generated 
with a Gaussian probability distribution that is uncorrelated 
from pixel to pixel and between different years. The noise 
variance per pixel is inversely proportional to the hit count 
for each individual year. This is expected to describe reason- 
ably accurately the expected noise behaviour of the actual 
WMAP data. The noise level impacts the exact value of the 
NILC coefficients (in particular at high galactic latitude and 
on small scales), and the exact size of the error bars in the 
measured CMB power spectrum. Small uncertainties in the 
noise, however, do not bias the estimated Ci. 



^ |http://www.apc.univ-paris7.fr/~delabrou/psm.htmI| 



Biases are investigated in the following way. Fifty dif- 
ferent realisations of the full simulations are generated, and 
processed in the exact same way as the actual WMAP data 
sets. The average systematic difference (i.e. bias) between 
the recovered power spectra and the input CMB power spec- 
trum is computed. The bias estimated in this way (displayed 
with green right-facing triangles in figure [BJ is typically less 
than 3 — 5% of the estimated Ci, and also typically smaller 
than the total errors on individual Ci bins (figure This 
average bias is nonetheless corrected for in our estimate of 
the CMB power spectrum (see sect ion IT4l for details). 

The observed bias is due to the combination of two ef- 
fects: additive residual foreground emission, which gener- 
ates a positive bias, and empirical correlations between con- 
taminants (foregrounds and noise) and CMB maps which 
gene rate the n egative ILC bias discussed at length in 
DelabrouiUeet al. (2009). 

Our de-biasing procedure deserves a discussion, since 
using this average bias estimated on simulations to de-bias 
(by subtraction) the actual CMB power spectrum obtained 
on real WMAP data is valid only if the simulations are rep- 
resentative of the real sky, i.e. 

• the CMB power spectrum assumed in the simulations 
must be close to the real one; 

• the simulated foregrounds must be representative of the 
sky (similar levels and behaviour) 

• the description of the instrument must be correct. 

The real CMB power spectrum and the real foregrounds are 
not known to infinite accuracy (otherwise there would be 
no point in trying to measure them again). However, the 
following points give us confidence in the representativeness 
of our simulations: 

• Concerning the CMB, a small error AC; in the simu- 
lated CMB Ci would generate a miscalculation of the multi- 
plicative ILC bias. The error in this bias estimate, however, 
is of second order (of order 2% times AC;), with a minor 
impact on the C; estimate. 

• The galactic foreground model in the PSM is built 
on the basis of the WMA P foreground analysis of 
iMiviUe-Deschenes et"al] (|2008| ). Even if the model has 
known limitations and uncertainties, it yields sky emission 
in the WMAP channels that is in good agreement with the 
observations. For safety in our power spectrum estimation 
however, a conservative mask being applied to exclude the 
galactic plane region from the calculation of temperature C; . 

• PSM point sources are based on real observations 
of radio sources with GB6 (ICawiser et all ligOSl). PMN 
llCriffith et al.ll 19941 : IWright et al.ll 19941: iGriflith et al.lll995l : 
IWright et al.lll996l). S UMSS (iMauch et al.ll2003l ) ar^NVSS 
( Condon et al. I 1 19981 ). and includes the full catalog of 
WMAP detected point sources. The PSM point source 
model is hence thought to be reasonably representative of 
the expected point source contamination in WMAP data 
sets. 

• The thermal and kinetic SZ effects are very faint as 
compared to WMAP sensitivity, and do not play a siginfi- 
cant role in this analysis. 

• Uncertainties connected to the description of the in- 
strument that are susceptible to impact our simulations 
are mostly calibration uncertainties and beam uncertainties. 



A needlet ILC analysis of WMAP 7-year temperature data 7 



VVWAP7yr i 
NILC i 



200 400 600 SOD 1000 1200 

Multipole{l) 




400 600 

Multipole(l) 



800 1000 




0.0000 




Figure 4. The red filled circles show the angular power spectrum 
estimated, after foreground subtraction with the NILC, using 7 
years of observations of WMAP. The blue filled squares show the 
7-year angular power spectrum published by the WMAP collab- 
oration. The black solid line shows the theoretical angular power 
spectrum for WMAP best-fit Lambda-CDM model. The top panel 
uses a linear scale in the horizontal axis, and the bottom panel a 
logarithmic scale. Our estimate is consistent with that of WMAP, 
although it has significantly less power at low multipoles [l < 15). 



Figure 5. The red filled circles show the difference between the 
angular power spectrum estimated using NILC maps and the the- 
oretical WMAP best fit angular power spectrum. The red open 
circles show the total error in the estimate of the angular power 
spectrum using our method. Similarly, the blue filled squares show 
the difference between the angular power spectrum published by 
the WMAP collaboration and theoretical best fit. The blue open 
squares show the error in the estimate of angular power spec- 
trum provided by WMAP. The top panel uses a linear scale in 
the horizontal axis, and the bottom panel a logarithmic scale. 



The impact of these, which turns out to be negligible as es- 
timated on simulations, will be discussed later on. 



4.4 Power spectrum estimation 

Residual noise in our CMB maps is inhomogeneous, primar- 
ily because of non-uniform sky coverage, with higher number 
of observations in the directions of ecliptic poles and rings 
around them of 45° radius. Taking into account this inhomo- 
geneity for computing the CMB power spectrum is expected 
to yield more precise estimates of the CMB Ci in multipoles 
where the noise is the main source of error. 

We follow the ge neral idea o f the n oise- weighting using 
needlets described in iFav et all (|2008l ). with a number of 
changes to adapt the method to our present analysis. Con- 
sider a bin in I, defined by a window Hi,{l) — [hi,{l)]^ re- 
stricted to the range in I limited by l'^ 



1(1+1) 



Hb{l) = 



2n (21 + 1) (Ai)b 





otherwise 



with ( AZ)(, = - -I- 1. The power Pj, of the CMB map 
filtered with hb{l) is an estimate of the CMB angular power 
spectrum Ce in the window Hb{l), with: 



1(1 + 1) 



(A/), 



2tv 



Ci 



(19) 



Now starting from a pair of reconstructed CMB maps for 
two different years of observation i and j, filtering each of 
them by hbil), we can form quadratic estimators of the band- 
averaged power spectrum between Zmin and /max: 

47r 



iVpix 



(20) 



where T(,,i(^) is the NILC estimate of CMB map for year i, 
filtered in the band hb{l), and A^pix is the number of pixels 
in the filtered map. 

In fact, for each single pixel ^, the quantity 



Pb'{0 = ^^TbAOTbAO 



(21) 
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is itself an estimator of Pt,. We form a noise- weighted aver- 
age of the estimators for all pixels and all pairs of maps of 
the form: 



iVpix A^yoar (iVy, 



1) 



using weights proportional to 
1 



47r 



Pb + Un 



(23) 



where A'yoar is the number of years of observation and 
o"noiso(C) is the variance of the pixel noise of the recon- 
structed 7-year NILC CMB map fihered by hb{l). The 
weights are normalised so that no power is lost: 



iVpix iVyoar (iVy, 



(24) 



5 i<j 



These weights a re very similar to those in equation 18 of 
I Fay et al] (|2008l ). used for for power spectrum estimation 
using a needlet frame, except that here we have a correction 
term o-'^omdV) / {^y^a^i ~ 1)5 coming from the fact that only 
cross spectra of the A^yoar maps are used in our estimator. 
When a large number of independent maps are used, this 
term is small and can be neglected (in our case, with A'yoar = 
7, the correction to the weights is of order 15% only). 

The value of cr^oisAO is estimated using 100 indepen- 
dent realisations of the reconstruction of the CMB by the 
NILC, and Pb is estimated by plugging in the WMAP best 
fit Ci in eauation ll9l Note that some imprecision on Pi, and 
fnoiso hS'S little impact on the estimation. In case of errors 
in the estimation of Pb and (jj^o;^^, the weights used are not 
perfectly optimal, but this induces no bias in the estimator. 

An apodised galactic mask Mioo^io" is used in the power 
spectrum estimation (see equation (|19|l and figure [Sjl. Such 
masking is equivalent to lowering the weight of some of the 
pixels (down to vanishing weights at the lowest galactic lat- 
itudes) . 

Figure |4] shows the binned NILC estimate of the angu- 
lar power spectrum of CMB temperature anisotropies, after 
subtracting the binned ILC bias. The estimate of the same 
obtained by the WMAP collaboration is plotted on the same 
figure for comparison. Figure [5] shows, for each of them, the 
difference with the theoretical angular power spectrum, and 
figure [6] displays the estimated NILC and calibration biases, 
compared to the statistical error of the estimator. 

The error bars in the estimated power spectrum include 
the statistical error of the estimator (noise and cosmic vari- 
ance terms) and the statistical error on the ILC bias term 
- the value of which is known only up to statistical error, 
by reason of the intrinsic variance of the empirical correla- 
tion between the CMB and the contaminants (foregrounds 
and noise). The power spectrum obtained using our analysis 
agrees well with that provided by the WMAP collaboration 
and is consistent with the ACDM model for WMAP best-fit 
cosmological parameters, except for a small but significant 
lack of power at low multipoles, as can be seen in figure 21 
We do not observe such systematic differences on our anal- 
ysis of simulated data sets. 
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Estimated bras (0.5% Calibration 
■ ■ Estimoted bios (ILC + Eoregroun 






600 
Multipale(l) 



Estimated bias (0.5% Calibration error) 
.. Estimoted bios (ILC + Foreground + Apodized mosk) 




Figure 6. The green right-facing triangles show the bias esti- 
mated by implementing the NILC on 50 WMAP-Iike simulated 
data. The yellow downward triangles show the estimate of the 
bias due to calibration error of the order of of 0.5%. The red 
open circles show the total error in the estimate of angular power 
spectrum using NILC. The top panel uses a linear scale in the 
horizontal axis, and the bottom panel a logarithmic one. 



4.5 Impact of calibration errors 

Calibration errors are a serious issue for precise separation 
of the C MB from foregroun ds using any type of ILC. As dis- 
cussed bv lDick et all (|2010l l. they ' conspire' with the ILC fil- 
ter to cancel out the CMB. This effect is particularly strong 
in the high signal to noise ratio regime, which is the case 
in our present analysis, in particular on large scales. We 
investigate the impact of this calibration bias by redoing 
the analysis using slightly modified calibration coefficients 
(differences of the order of 0.5%), and computing the dif- 
ference between the CMB spectra estimated in both cases. 
The result, plotted in figure [S] shows that here even at large 
angular scale, the bias due to calibration error is small. We 
hence neglect this effect in our estimate of errors on Ce - Note 
however that this calibration bias will be a concern for those 
experiments which measure CMB with higher signal to noise 
ratio than WMAP. 



4.6 Impact of beam uncertainties 

The uncertainty on beam shapes is equivalent to a calibra- 
tion error which depends on the harmonic mode I. As this 
uncertainty is essentially at high /, where the signal to noise 
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o Error ottribuied to cosmic vorionce 
A Error ottributed to measurement errors 

□ Error attributed to bias(lLC + Foreground + Apodized mask) 
o NILC error 



Table 2. Comparison of NILC estimate of binned angular power 
spectrum (after bias correction) with the estimate of binned an- 
gular power spectrum provided by WMAP collaboration. 
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Figure 7. The red circles sliow the total error in the estimate 
of the CMB angular power spectrum using the NILC. The black 
open diamonds show the error due to cosmic variance. The blue 
upward triangles show the measurement error. The sky blue open 
squares show the error due to bias uncertainty. The horizontal 
axis of the top plot is scaled linearly, and that of the bottom plot 
is scaled logarithmically. 



ratio is worse than on large scales, the bias, for small beam 
shape errors, is not expected to impact much the CMB re- 
construction here, and is also neglected in the present anal- 
ysis. 



5 DISCUSSION 

5.1 Low power on large scales? 

The temperature power spectrum we obtain seems to be, on 
large scales, systematically lower than the theoretical 'best 
fit' power spectrum (for £ < 15) and the WMAP measure- 
ment (for i < 40, except for the quadrupole, see figure IH 
bottom left panel, and table [3]). 

Considering that error bars on these scales are cosmic 
variance dominated, the difference between the power spec- 
trum measured by WMAP and by us from the same original 
data set is problematic, and deserves a discussion and ad- 
ditional investigation. A systematic shift in the estimated 
power spectrum at all I below 40 is bound to impact the 
interpretation of the observations, and the inferred values 
and limits for cosmological parameters. In the following, we 
investigate various possible origins for the discrepancy. 
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5.1.1 Bias 



to the method of spectrum estimation? 



Our method is based on the computation of the needlet 
transform of a masked sky. We check that the method itself 
is not biased, by computing in the exact same way the power 
spectrum of 100 pure CMB maps, randomly generated us- 
ing different seeds. The average estimated power spectrum 
is computed, and individual recovered power spectra for the 
100 realisations visually inspected. We find no evidence of a 
bias in our power spectrum estimate. 
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Table 3. Difference between the NILC estimate of binned angular 
power spectrum and the angular power spectrum provided by the 
WMAP collaboration. Except for the quadrupole, for which our 
estimate is larger (after debiasing) and except the bins centred at 
I = 9 and I = 27, our estimated power spectrum is lower than the 
WMAP one. 
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Figure 8. The red filled circles show the average estimated an- 
gular power spectrum using NILC from WMAP 7-year-like sim- 
ulations. The black solid line show the theoretical angular power 
spectrum for WMAP best-fit ACDM model. The top panel uses 
a linear scale in the horizontal axis, and the bottom panel a log- 
arithmic one. 



5.1.2 Bias due to the complete analysis chain? 

At a next level, we check whether our complete data process- 
ing method can generate a bias. We generate two times 50 
simulations of the NILC maps for the 7 independent years. 
The first 50 coincide with the simulations used to estimate 
an average ILC bias (green curve in figure |6} and correct for 
it in the real data (the effect is small). The other 50 sim- 
ulations are used as test data, for which we implement our 
complete analysis pipeline exactly as done on the WMAP 
data. For each of them, in particular, we de-bias the result- 
ing power spectrum using the average bias in the first 50 
simulations, in the exact same way as was done on the real 
data. 

Figure[8]shows the average angular power spectrum (af- 
ter bias correction) of CMB temperature anisotropics on 
these test data sets. The power spectrum matches extremely 
well the theoretical input angular power spectrum, and does 
not show any lack of power at large angular scale. This 
demonstrates that our method does not generate biases by 
itself. Nor the needlet ILC, nor the way we implement our 
power spectrum estimation can be responsible for a bias - 
assuming our data model is right. 

We note that the bias correction at / = 2 is of more 
than 3 x 10~* mK'^, and that the error on the estimate of 
this bias is itself of order 5 x 10"'' mK'^. This can explain 
the difference between the WMAP measured quadrupole, 
and our estimate. 

5.1.3 Are simulations representative? 

The representativeness of the simulated data sets is of course 
a major concern. For instance, if our simulated maps had 
exceedingly large galactic foregrounds, the average bias on 
the simulations due to residual foregrounds would be too 
large, and subtracting this residual from the real data would 
result in an over-correction, and hence a negative bias. 

Although significant effort has been put to generate sky 
simulations as realistic as possible, such a bias due to mod- 
elling errors cannot be fully ruled-out (the modelling uncer- 
tainty is not easily estimated). For this reason, we look for 
confirmation using other CMB maps and different analyses. 

5.1.4 Confirmation using other WMAP CMB maps? 

We have also estimated the power spectrum of WMAP maps 
directly using the original V and W frequency channel maps, 
the foreground-reduced V and W channel maps published by 
the WMAP collaboration, and the WMAP ILC map itself, 
applying different galactic masks. These estimates also result 
in low £ CMB multipoles. In figure [§] we compare the mul- 
tipoles computed directly on the 7-year NILC map and on 
the 7- year WMAP ILC. They are almost indistinguishable 
at low I. The very low value of the spectrum of the difference 
between those two maps shows that the difference between 
our present power spectrum and the published WMAP spec- 
trum is not due to significant differences in the large scale 
modes of the maps themselves, but on the power spectrum 
estimation method (including debiasing). 

Power spectra computed directly from the 7- year V and 
W channels also tend to be lower than the WMAP theoret- 
ical best fit. A fraction of the deficit in power on the largest 
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scales (/ = 2 to / = 4) is due to the mask (there is no bias 
correction here, except for a global /sky factor). In most of 
the measured multipoles for I < 30, there seems to be less 
power in the maps than theoretically predicted. This is par- 
ticularly true for conservative masks (M2o° ,20° and kp2). We 
note that the power in the V band map decreases with in- 
creasing galactic cut. For the M2o°,20° mask, the spectra of 
all maps (Needlet ILC, WMAP ILC, V band and W band) 
are almost indistinguishable, and on average significantly 
lower than the theoretical best fit, for 5 < I < 25. 



5.2 Comparison with other work 

We have not been able to identify a source of bias in our anal- 
ysis which would explain this difference between our power 
spectrum and the one published by the WMAP collabora- 
tion. A significant amount of connected work has been done 
by various authors to extract a CMB map, and compute 
CMB power spectra from WMAP data, with different tech- 
niques. We now discuss how our present work differs from 
such previous analyses with similar objectives, and compare 
the results obtained with various methods and on various 
W MAP data releases ( i n part icular on large scales). 

iDelabrouille et all \200^ make a CMB map from 
WMAP 5 year data on needlet (spherical wavelet) domains. 
Our present work is an extension implemented on WMAP 7- 
year data. As an additional analysis step, we also use the cor- 
responding NILC weights to obtain clean CMB maps from 7 
individual years, and compute cross spectra to estimate the 
CMB power spectrum. 

A CMB power spectrum has bee n published with each 
of the various WMAP data releases ( Hinshaw et al.l [20031 . 



l200'7l : lNolta et al.ll2009l : lLarson et al.ll201ll ). The CMB power 
spectrum on large scale is unexpectedly low in the anal- 
ysis of the first year data. In subsequent releases, only the 
quadrupole remains significantly lower than the WMAP best 
fit model spectrum. 

The most recent temperature power spectrum produced 
by the WMAP collaboration is obtained on the basis of two 
different techniques for large and small scales. For £ ^ 32, 
the spectrum is obtained using a Blackwell-Rao estimator 
applied to a chain of Gibbs samples based on the seven-year 
ILC map masked with the KQ85y7 mask, which cuts 21.7 % 
of the sky close to the galactic plane. The difference in power 
spectrum at low multipoles between this estimate and ours 
must originate from either the exact way the power spectrum 
is estimated on large scales (including debiasing) , but we do 
not know at present whether the observed difference is to be 
expected on the basis of differences in the analysis pipeline, 
nor its exact significance, as this requires checking details 
of both analyses, some of which are not available to us at 
present. 

Among the differences, we note that our sky model is 
based on 4 galactic components (synchrotron, free-free, spin- 
ning dust and thermal dust), with spectral indices vary- 
ing over the sky for synchrotron and thermal dust, while 
the WMAP 3-year analysis, for instance, uses a three- 
component model of sky emission. It is possible that a 
bias correction based on simulations with a 3-component 
ISM emission yields a higher Ci estimate than with a 4- 
component one. If the origin of the discrepancy is there, it 
emphasises the need for complete models of sky emission, in 
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Figure 9. CMB power spectra obtained from different WMAP 
maps. The red circles show the angular power spectrum obtained 
by direct calculation of tlie spectrum of the NILC CMB obtained 
on WMAP-7year band average maps, masked with an apodised 
Galactic mask. The mask used is Mjqo jgo (same mask as used 
in our analysis) for the top panel, M200 20° foi" the middle panel, 
and kp2 (including point sources) for the bottom panel. The 7- 
year NILC map power spectrum is fully consistent with the power 
spectrum of the WMAP 7- year ILC map (blue squares), masked 
in the same way. No correction for the ILC bias (which depends 
on the mask used) was applied here. The yellow diamonds show 
the angular power spectrum of the difference of the two ILC maps. 
This difference, very low on large scales, rises at higher / because 
of the larger noise residuals in the WMAP ILC map. The down- 
ward and upward triangles show the power spectra of the 7-year 

V and W channels respectively. As can be seen by comparing the 
various panels, galactic foregrounds contaminate significantly the 

V channel at galactic latitudes between 10 and 20 degrees (and, 
to a lesser extent, the W channel as well). Note also that the 
lowest multipoles depend somewhat on the exact mask used, as 
expected for 'pseudo' spectra. 
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which the complexity of the sky can be varied at the level 
of present modelling uncertainties. 

Another particularity of our analysis is the use of cross 
spectra between data from different years of observations, 
which in principle makes our analysis impervious to large 
scale residuals of low frequency noise correlated between 
the various WMAP channels. Such hypothetical residuals, 
however, are unlikely to be the main reason for the low-/ 
discrepancy, as demonstrated by the comparison of the 7 
year NILC map and the WMAP ILC map at low L 

At higher multipoles, the WMAP spectrum is estimated 
on the seven-year, template-cleaned V- and W-band maps, 
using (pseudo) cross-spectra between maps obtained for in- 
dividual differencing assemblies and individual years. Above 
£ = 600, the maps are inverse-noise weighted to reduce the 
noise variance in the final power spectrum estimate. The 
WMAP power spectrum is in good agreement with our es- 
timate. 

In addition to the analyses published by the WMAP col- 
laboration, the CMB power spectrum inferred from WMAP 
data se ts has been investigated by a number of authors. 

.Eriksen et all (|2007l 'l. in particular, partly re-analyse 
WMAP 3-year data sets with different methods. They find 
large scale power spectra systematically lower than the 
WMAP best fit and published power spectrum, although not 
as low as our present result for 3 ^ Z ^ 6. All the estimated 
spectra they publish are somewhat different, a feature they 
interpret, on the basis of differences when the galactic mask 
is extended, as plausibly due to contamination by residual 
galactic foreground emission. 

On the foreground cleaning side, CMB reconstruction 
using ILC method s has been performed by several authors. 
ISaha et all l|2008l ) have implemented an ILC on WMAP 1- 
year and 3-year data in the harmonic domain. The primary 
goal of their work is to compute the CMB power spectrum, 
rather than producing a clean CMB map. They also analyse 
the bias produced at low / by the ILC in the measured power 
spectrum. This bias is estimated in the case of an imple- 
mentation in the harmonic domain, and an equivalent bias 
correction is made in our cas e on the basis of Monte-Carlo 
simulations. In their analysis, [Saha et al.1 Hooi) also find, as 
we do in our present analysis, that both the quadrupole and 
the octupole are significantly lower than the WMAP best 
fit. 

ISamal et al.1 (|2010h estimate CMB polarisation and 
temperature power spectra using linear combination of 
WMAP 5 year maps. The authors use different combina- 
tions of individual WMAP differencing assemblies, rather 
than different years, to compute independent sky maps, but 
find, again, somewhat lower CMB power on large scale than 
both the WMAP estimate and best fit. 

Hence, significant differences between the low CMB 
multipoles are found by different authors, using different 
methods. All the analyses are based on the same sky emis- 
sion, and hence should agree to up to errors which exclude 
cosmic variance. Even worse, these analyses sometimes start 
from the same original data, so that even the instrumental 
noise errors should be the same. This demonstrates that 
the low I power spectrum estimation is sensitive to details 
of the analysis pipeline. The exact origin of the discrep- 
ancy between the various estimates remains to be under- 
stood, but the fact that there still is debate on the exact 



value of the low I CMB multipoles is a concern, which il- 
lustrates the difficulty of component separation and power 
spectrum estimation on large scales. Solving these issues is 
important for the tentative measurement of primordial ten- 
sor modes with Planck and with future CMB polarisation 
experiments such as the r ecently proposed CMBPol, EPIC, 
and COrE space missions (iBaumann et aL l2009l : iBock et al.l 
l2008l : lThe COrE Collaboration et al.ll201ll ). 



6 CONCLUSION 

The precise measurement of cosmological parameters has set 
a new goal in the effort to estimate the angular power spec- 
trum of the CMB by present and future CMB experiments. 
This effort is primarily boosted by the ever increasing im- 
provement in the sensitivity and resolution of the CMB ob- 
servations. However, the removal of foreground contamina- 
tion from the observed sky maps is a mandatory preliminary 
step for the accurate estimation of angular power spectrum. 

In this paper, we have described a new methodology to 
estimate the angular power spectrum of CMB temperature 
anisotropy from WMAP 7-year data. We have used linear 
combination of sky maps decomposed on a frame of spherical 
wavelets (needlets), to construct a map of CMB anisotropics 
with low contamination from foreground signals and instru- 
mental noise. The CMB temperature anisotropy map has 
been estimated by implementing the NILC on WMAP 7- 
year band average maps at the highest resolution (with 
HEALPix pixelisation parameter nside equal to 1024). We 
have included, in our analysis, three foreground templates 
in the set of analysed observations for better performance 
of our component separation. Our method for CMB clean- 
ing does not rely strongly upon any assumed model of fore- 
ground emission and detector noise properties, and hence, 
is not very sensitive to the uncertainty and insufficiency in 
foreground and noise modelling. However, the cleaned CMB 
map always contains, at some level, non-vanishing residual 
foreground and residual noise, superimposed on actual CMB 
signal. 

To minimise the impact of uncertainties in detector 
noise levels and correlation between channels, the angular 
power spectrum of the CMB has been estimated from all 
possible (twenty one in our case) cross power spectra of 
clean CMB maps for 7 individual years of observations with 
independent noise. The estimates of CMB maps for the 7 
individual years are obtained by using the NILC weights 
obtained by implementing a needlet ILC on WMAP 7- year 
band average maps. 

To reduce the impact of foreground signals on the mea- 
sured CMB power spectrum, in addition to using low fore- 
ground linear combinations of input WMAP channels and 
of ancillary data, we have applied the point source mask 
provided by the WMAP collaboration, and then filled the 
masked regions by a interpolation procedure. We have also 
used an apodised symmetric mask to lower the impact of 
residuals emission from the interstellar medium in our anal- 
ysis. 

Biases due to random correlation of the CMB with the 
ISM in our NILC pipeline have been estimated using 50 re- 
alisations of WMAP-like simulations. These biases turn out 
to be smaller than the other sources of error on most of the 
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range of I, but not completely negligible, so we use the av- 
erage bias measured on these 50 simulations to de-bias our 
estimate on WMAP data. We find that our de-biased angu- 
lar power spectrum agrees well with the estimate of angular 
power spectrum provided by the WMAP collaboration, ex- 
cept at large angular scale, where our estimate seems to be 
systematically lower, but where the fundamental uncertain- 
ties are high. A number of tests, performed on simulated 
data sets, using different approaches to power spectrum es- 
timation on WMAP data, confirm that the released WMAP 
maps seems to contain somewhat less CMB anisotropics on 
large scales than expected from the WMAP best fit model 
(marginally within cosmic variance errors, as can be seen in 
the bottom panel of figure[5l). This still has to be elucidated. 

The error bars in our estimate of de-biased angular 
power spectrum have been computed very carefully. They 
comprise estimates of the measurement and cosmic vari- 
ance error on the basis of the internal scatter of 21 inde- 
pendent cross-spectra, and of the error in the estimate of 
the ILC bias. Errors due to miscalibration or to imperfect 
beam knowledge have been shown to be small enough to 
be negligible, considering the calibration uncertainty esti- 
mates provided by the WMAP collaboration. Our total error 
bars are comparable to the error in the estimate of angular 
power spectrum provided by the WMAP collaboration. As 
the measurement errors are estimated from the internal scat- 
ter of independent cross spectra for each bin in I, they do 
not rely on a model of WMAP detector noise. 

The fact that there still is not a convincing consensus 
in the scientific community on the exact value of the low / 
CMB multipoles is a concern, which demonstrates the diffi- 
culty of extracting precisely the CMB emission from galactic 
foregrounds on large scales. 
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